1.85 Approximation for Min-Power Strong Connectivity 



G. Calinescu * 



in 



q 

o 



(N 

Abstract 

o ; 

fSj , Given a directed simple graph G = (V,E) and a cost function c : E — )• R + , the power 

of a vertex w in a directed spanning subgraph H is given by ph{u) = m a x rae-E(#) c(w), and 
corresponds to the energy consumption required for wireless node u to transmit to all nodes v with 
uv G E(H). The power of if is given by p(ii) = ]Cuev Ph{u). 

Power Assignment seeks to minimize p(H) while if satisfies some connectivity constraint. 
In this paper, we assume E is bidirected (for every directed edge e £ E, the opposite edge 
exists and has the same cost), while if is required to be strongly connected. This is the original 
power assignment problem introduced by Chen and Huang in 1989, who proved that a bidirected 
minimum spanning tree has approximation ratio at most 2 (this is tight). In Approx 2010, we 
introduced a Greedy approximation algorithm and claimed a ratio of 1.992. Here we improve 
the analysis to 1.85, combining techniques from Robins-Zelikovsky (2000) for Steiner Tree, and 
Caragiannis, Flammini, and Moscardelli (2007) for the broadcast version of Power Assignment, 
together with a simple idea inspired by Byrka, Grandoni, RothvoB, and Sanita (2010). 
| The proof also shows that a natural linear programming relaxation, introduced by us in 2012, 

' has the same 1.85 integrality gap. 

ON ' 
CO 

CO ■ 

^ 1 Introduction 

o 

CN ; There has been a surge of research in Power Assignment problems since 2000 (among the earlier 
T"! ■ papers are lfT8l 1221 [T3l ) This class of problems take as input a directed simple graph G — (V,E) and 
. ^ . a cost function c : E — )■ R + . The power of a vertex u in a directed spanning simple subgraph H of 
^ \ G is given by Ph(u) = max uv( z E ( H ) c(uv), and corresponds to the energy consumption required for 
wireless node u to transmit to all nodes v with uv E E(H). The power (or total power) of H is given 
by p(H) = J2ugvPh(u)- 

The study of the min-power power assignment was started by Chen and Huang (S), which consider, 
as we do, the case when E is bidirected, (that is, uv E E if and only if vu E E, and if weighted, 
the two edge have the same cost; this case was sometimes called "symmetric" or "undirected" in the 
literature) while H is required to be strongly connected. We call this problem Min-Power Strong 
Connectivity. We use with the same name both the (bi)directed and the undirected version of G. (SI 
prove that the bidirected version of a minimum (cost) spanning tree (MST) of the input graph G has 
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power at most twice the optimum, and therefore the MST algorithm has approximation ratio at most 
2. This is known to be tight (see Section©. 

We improve this to 1.85 by combining techniques from Robins-Zelikovsky [fT9l l20l for Steiner 
Tree, Caragiannis, Flammini, and Moscardelli [|6) for the broadcast version of symmetric Power As- 
signment (assuming a bidirected G = (V, E, c) and a "root" u G V is given, H must contain a directed 
path from u to every vertex of G), together with a simple idea inspired by Byrka, Grandoni, RothvoB, 
and Sanita (2010). 

Very restricted versions of Min-Power Strong Connectivity have been proven NP-Hard lTT5l ITOl mi . 
Other than £4), we are not aware of better than a factor of 2 approximation except for [7], (where 
c : E — > {A, B}, for < A < B; see also flU), [|2l (where c is assumed to be a metric), and the 
exact (dynamic programming) algorithms of [fl5l for the specific case where each vertex of G maps 
to a point on a line, and c(uv ) is an increasing function of the Euclidean distance between the images 
of u and v. A related version, also NP-Hard, asks for H to be bidirected (also called "undirected" 
or "symmetric" in previous papers). This problem is called Min-Power Symmetric Connectivity, and 
the best known ratio of 5/3 + e HI is obtained with techniques first applied to Steiner Tree; when 
c : E — )■ {A, B} one gets 3/2 with the same method [17J. In fact, many but not all power assignment 
algorithms use techniques from Steiner Tree variants (or direct reduction to Steiner Tree variants; 
these connections to Steiner Tree are not obvious and cannot be easily explained), and in particular 
Caragiannis et al [6] uses the relative greedy heuristic of Zelikovsky [23 ] . New interesting techniques 
were also developed for power assignment problems, as in fH6ll . an improvement over [[141. 

The existing lower bound of the optimum, which we use, is the cost of the minimum spanning tree 
of G. Indeed (argument from [8]), the optimum solution OPT contains an in-arborescence rooted at 
v, for some v E V, and then, for all u G V \ {v}, Popt(u) is at least the cost of the directed edge 
connecting u to its parent in this in-arborescence, whose total cost is at least the cost of the minimum 
spanning tree of G. 

We also use a relative greedy method as in [1231 [191 ; Robins-Zelikovsky lfT9l is rarely used as a 
technique, and not by only citing the ratio (improved by now in Oj). We use the natural structures 
of j6l to improve over the minimum spanning tree: these are stars, directed trees of height 1. Our 
second lower bound (improved over [4]) comes from "covering" the edges of the a spanning tree by 
the stars of the optimum solution. With precise definitions later, we just mention that an edge of a 
tree is covered by a star if it is on a path of the tree between two vertices of the star. A "cheap" 
fractional covering can be easily obtained from either optimum or the linear programming relaxation. 
In our earlier work 0], we used an integer cover which was extremely hard to obtain. Using fractional 
covers (inspired by [0) is the only significant difference of this version versus [4J. Also, interestingly, 
the submodularity of the covering function is only used implictly. 

2 Preliminaries 

In directed graphs, we use arc to denote a directed edge. In a directed graph K, an incoming arbores- 
cence rooted at x G V(K) is a spanning subgraph T of K such that the underlying undirected graph 
of T is a tree and every vertex of T other than x has exactly one outgoing arc in T. 

Given an arc xy, its undirected version is the undirected edge with endpoints x and y. Arcs xy 
and yx are antiparallel, and the antiparallel arcs resulting from undirected edge uv are uv and vu; if 
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Figure 1 : Tight example for the performance ratio of the MST algorithm. In both cases, the solution 
is bidirected and the undirected version of the arcs of the solution is given by solid edges. For each 
vertex, its power in the solution is written next to it. (a) The MST-based power assignment needs total 
power 2n. (b) Optimum power assignment has total power n(l + e) 2 + (n — l)e 2 + 1 — Y n + 1, 



undirected edge uv has cost, then each of the two antiparallel arcs resulting from undirected edge uv 
have this cost. We sometimes identify a spanning tree T with its set of edges. 

An alternative definition of our problem (how it was originally posed) is: we are given a simple 
undirected graph G = (V, E) and a cost function c : E — >■ R + . A power assignment is a function 
p : V — > R + , and it induces a simple directed graph H(p) on vertex set V given by xy being an arc of 
H(p) if and only if {x, y} G E and p(x) > c({x, y}). The problem is to minimize Ylu&v P( u ) subject 
to H(p) being strongly connected. To see the equivalence of the definition, given directed spanning 
subgraph H, define for each u E V the power assignment p(u) = Ph{u). 

The following known example (see Figure [TJ) shows that the ratio of 2 for the MST algorithm is 
tight. Consider 2n points located on a single line such that the distance between consecutive points 
alternates between 1 and e < 1, and let the cost function c be the square of the Euclidean distance 
Then the minimum spanning tree MST connects consecutive neighbors and has power p(MST) = 2n. 
On the other hand, the bidirected tree T' with arcs connecting each other node (see Figure Eb)) has 
power equal p(T') = n(l + e) 2 + (n — l)e 2 + 1. When n — > oo and e — > 0, we obtain that 
p(MST) /p(T') — > 2. On the other hand (argument taken from [|8]), for any input graph, the power of 
the bidirected minimum spanning tree T is at most 

p(T) = > max c(vu) < > > c(vu) = 2c(T) < 2 opt 

^— ' u\vu&E(T) > - S-^ '- 

veV v&V u\vu&E{T) 

where opt = p( OPT) for an optimum solution OPT (the last inequality is from the introduction). 

The example above may give intuition on how Power Assignment (and even more specifically, 
Min-Power Symmetric Connectivity - the variant when H must be bidirected mentioned in the in- 
troduction) relates to the ^-restricted Steiner trees, with stars (trees of height 1) taking the place 
of restricted components. Another example from [fl], (see Figure |2] and the following paragraph), 
shows how Min-Power Strong Connectivity differs from Min-Power Symmetric Connectivity, and 
may give intuition how Min-Power Strong Connectivity relates to Travelling Salesman, and also Min- 
Cost Strong Connectivity and Min-Cost Two-Edge Connectivity (a two-edge-connected graph has an 
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Figure 2: Total power for Min-Power Strong Connectivity can be half the total power for Min-Power 
Symmetric Connectivity. In both cases, the solution is represented by solid segments (arrows meaning 
not bidirected), and next to each vertex its power is given, (a) Minimum power assignment ensuring 
strong connectivity has total power n + n 2 e 2 = n + n 2 ^ = n + 1. (b) Minimum power assignment 
ensuring symmetric connectivity has total power (2n — 2) + [n{n — 1) + 2)e 2 = 2n — 1 — - + the 
solution is bidirected and the undirected version of the arcs of the solution is given by solid segments. 

edge orientation that makes it strongly connected - see for example Chapter 2, written by A. Frank, of 
ifTTlO . However we cannot think of direct reductions either way, and, as we mention in Conclusions, 
the methods we use only apply to certain instances of Min-Cost Strong Connectivity and Min-Cost 
Two-Edge Connectivity. 

The power of a Min-Power Strong Connectivity optimum solution can be almost half the power 
of a Min-Power Symmetric Connectivity optimum solution for the same instance: we present a series 
of examples illustrated in Figure |2] The n(n + 1) vertices are embedded in the plane in n groups of 
n+1 points each. Each group has two "terminals" (represented as thick circles in Figure |2]), and the 
2n terminals are the corners of a regular 2n-gon with sides of length 1 . Each group has another n — 1 
equally spaced points (dashes in Figure |2]) on the line segment between the two terminals. The cost 
function c is the square of the Euclidean distance. It is easy to see that a minimum power assignment 
ensuring strong connectivity assigns a power of 1 to one thick terminal in each group and a power of 
e 2 = (l/n) 2 to all other points in the group - the arcs going clockwise. The total power then equals 
n + 1. For symmetric connectivity it is necessary to assign power of 1 to all but two of the thick 
points, and of e 2 to the remaining points, which results in total power (2n — 2) + (n(n — 1) + 2)e 2 = 
2n — 1 — l/n + 2/n 2 . Also, keep in mind that the minimum spanning tree solution is a symmetric 
solution. 

3 The Approximation Algorithm 

This section is dedicated to proving the main result of this paper: 

Theorem 1 There exists a polynomial time algorithm for Min-Power Strong Connectivity with ap- 
proximation ratio 1.85. 

The outline of the proof of Theorem [TJ is as follows: 
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1. In the first subsection, we present the greedy algorithm, preceded by necessary notation. 

2. Then we establish (LemmaO that the algorithm's output is strongly connected. 

3. A second subsection gives the analysis. It starts with the new lower bound, the fractional cover 
of tree edges by stars. 

4. Finally, Lemma [5] shows how the Greedy method of [19] (applied with new parameters) com- 
bines the two lower bounds (the one above, and the cost of the minimum spanning tree) to 
obtain the claimed approximation ratio. 

3.1 The algorithm 

Our algorithm uses a greedy approach similar to [f23l H9l |6|. Let T be the undirected minimum 
spanning tree of G. The fact that T has minimum cost will be not further used, except to note opt > 
c(T), where opt = p( OPT) for an optimum solution OPT. Let T be the bidirected version of T. 

For u E V and r E {c(uv ) \ uv E E}, let S(u, r) be the directed star with center u containing all 
the arcs uv with c{uv ) < r; note that r is the power of S, also denoted by p(S). For a directed star S, 
let E(S) be its set of arcs and V(S) be its set of vertices. 

For given S(u, r), let Q(u, r) be the set of edges e of T such that there exist x,y E V(S(u, r)) 
with e on the path from x to y in T. Let Q(u, r) be the set of arcs e of T such that there exist 
x E V(S(u, r)) with e on the directed path from u to x in T; it is easy to verify that the undirected 
version of Q(u, r) is Q(u, r). 

For a collection A of directed stars S(ui, r^), define Q(A) = \Js( Ui n)eA Q( u ^ r «) an( ^ f(-A) = 
J2e£Q(A) c ( e )- ^ e sometimes write Q(S) instead of Q({S}. The function f(A) is known to be 
monotone and submodular (see an example in [2T), pages 768-769). For S = S(u,r), recall that 
fA(S) = f(A U {S}) - f(A) = J2eeQ(u,r)\ Q (A) c ( e ) = £ e e^(s) c ( e )> where ^(S) is defined to be 
those arcs of Q(u, r) for which the undirected version is not in Q(A). 

The algorithm starts with M = T as the set of arcs, and adds directed stars to collection A (initially 
empty) replacing arcs from M to greedily reduce the sum of costs of the arcs in M plus the sum of the 
powers of the stars in A. For intuition, we mention that this sum is our upper bound on the power of 
the algorithm's output. To simplify later proofs, the algorithm makes changes (adding directed stars 
and removing arcs from M) even if our sum stays the same. Assume below that 0/0 = 1. To be 
precise: 

Algorithm Greedy: 

A^Q,M^f 

While (f(A) < c(T)) do 

(u,r) <- argmax (u , r , ) /^(S'(w / ,r / ))/r / 

M <- M\I A (S(u,r)) 

A <- AU {S(u,r)} 
Output \J SeA E(S) U M 
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Each of figures |3]and|4] shows two iterations of the algorithm. For intuition, we mention that this 
algorithm "covers" undirected edges of the minimum spanning tree by "stars" and when implemented, 
it is a variant of Chvatal's [9] greedy algorithm for Set Cover. 




Figure 3: Costs are not relevant here, (a) Initial M. (b) Solid arcs give [j s&A E{S) U M after adding 
to A the star Si (with thick solid arcs) centered at x. The arcs removed from M are dashed, (c) 
Solid arcs give Us&A E(S) U M after adding to A the star S2 (with thick solid arcs) centered at y. 
The arcs removed from M are dashed. Note that the algorithm does not remove arc uv. Once an arc 
from a pair of antiparallel arcs of T is removed, the algorithm keeps the other arc, since there are 
cases (as in Figure H]) when not all arcs of Q(u, r) can be removed while keeping strong connectivity, 
and benefiting from removing arcs whose antiparallel arc has already been removed from T (when 
possible) destroys the "submodularity" implictly needed in the approximation ratio proof. 

Fact 2 Note that, unless f(A) = c(T), a star S(u,r) always exists for which f^(S(u, r)) > and 
fA(S(u, r))/r > 1. Indeed, as long as a pair of antiparallel arcs e' and e" are in M, we can pick as 
next star S(u, r) the one given by u being the tail of V and r = c(e'). 

Thus, as written, the algorithm can have iterations that do not change the output, i.e. above the 
star S(u, r) above could have just the edge e' and be added to A while e' is removed from M. 
Lemma 3 The output of Greedy is a spanning strongly connected subgraph ofG. 
Proof. We prove the following invariant: X := [J S( z A E(S) UM gives a spanning strongly connected 
subgraph whenever the while condition is checked by the algorithm. Moreover, suppose we remove 
from T all edges for which both antiparallel arcs appear in M, splitting T in components with vertex 
sets Tj, for some range of i. We prove that for every i and every x,y E Tj, there exists a directed path 
P from x to y using only vertices of Tj and arcs from X. 

Proving the invariant holds is done as always by induction on the number of iterations. The 
invariant is true before the first iteration, when each Tj has just one vertex, so consider the moment 
when a star S = S(u,r) is added to A. Figure |5] may provide intuition. We add arcs uz, for z E 
V(S) \ {u}, while removing from M (and from X) arcs xy if yx E M and there is some z such that 
xy is on the directed simple path from u to z in T. The same effect is obtained if we do this change 
for each z E V(S) \ {u} one after the other, instead of all such z at the same time. 

Let P be the simple path in T from u to z, and let Xj?/j, for 1 < i < k, be, in order, the arcs 
of M on P such that also yiXi E M. Thus the change to X consists of adding the arc uz and 
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(a) (b) (c) 

Figure 4: Costs are not relevant here, (a) Initial M. (b) \J SeA E(S) U M after adding to A the star Si 
(with thick solid arcs) centered at x, and removing arcs from M. (c) Solid arcs give {J SeA E(S) U M 
after adding to A the star S 2 (with thick solid arcs) centered at y. The arcs removed from M are 
dashed. Note that the algorithm cannot remove arc uv even though uv G Q(S 2 ), since u and v 
become disconnected. 

removing all arcs x^; note that if k = no arc is removed and our induction step is complete. Let 
M' = M \ {xiyi, . . . , x k y k } and X' — X \ {xiyi, . . . , x k y k } U {uz}. We need to show that X' and 
M' satisfy the conditions from the induction hypothesis. 

Let us split T into components by removing all the undirected edges xy with both antiparallel arcs 
xy and yx in M; in particular all Xiyi, for 1 < i < k, resulting in components Tj. 

By induction, X contains the following directed paths: Pi from xi to u, P 2 from x 2 to yi, 
P k from x k to yfc_i, and P k+ i from z to y k , and each of these paths stays in the same component Tj 
of the split of T by M - as in the previous paragraph. Thus none of these paths uses Xiyi or y t Xi. 
Putting together these paths, the arcs yiX { , for 1 < i < k, and the arc uz, we have a directed cycle C 
containing none of the arcs x^, for 1 < i < k. Any arc removed can be replaced, when discussing 
connectivity, with a path around the cycle C, and so (V, X') is strongly connected, as required. 

We now split T into components by removing all the undirected edges xy with both xy and yx 
in M', obtaining components T(. Note that none of Pi, 1 < % < k + 1, from above, has an arc with 
endpoints in two distinct components T[ (as T[ is the union of several Tj). As all the edges on the 
path from u to z in T do not have anymore both antiparallel arcs in M' , all the vertices on this path 
including u, z and all x i: yi axe in the same component T'- of the split of T by M' . Thus all the arcs of 
C have their endpoints in the same component of the split of T by M' . 

We prove that for every % and every x, y G T[, there exists a directed path P' from x to y using 
only vertices of T[ and arcs from X' . First, let us describe a path P from x to y using only arcs of X: 
find the path from x to y in T, and let ZiWi, for 1 < ? < q, be, in order, the arcs of M on P such that 
also G M. When g = 0, x, y G 1} for some j (same component of the split of T by M) and by 
induction, a path P from x to y exists in X using only vertices inside Tj. We pick P' = P, and indeed 
P' only uses arcs of X' since the arcs of X \ X' (same set as M \ M') cross from one component to 
another of the split of T by M. Assume now q > 0. Notice that all unordered pairs ZiWi belong in 
the set of unordered pairs xjyj on the simple path from u to v mentioned earlier, or else we cannot 
have that x and y belong to the same T( of the split of T by M'. Also, by induction, X contains the 
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Figure 5: Rounded rectangles show the components T it dashed before S is added to A (the split 
of T by M) and solid afterward (the split of T by M'). Arcs of M crossing from one component 
to another are given by thin arcs, S by the four thick arcs uzi ) uz 2) uz 3 , uz^, and the dashed arcs are 
those removed from M when S is added. 

following directed paths: P x from x to z\, P 2 from w 1 to z 2 , . . . , P q from u> 9 _i to z ? , and P 9+1 from 
w q to y, and each of these paths stays in the same component of the split of T by M. Thus none of 
these paths uses ZiWi or WiZi. Next, obtain P' by replacing in P, if necessary, arcs of X \ X' (same 
set as M \ M ') by arcs of C, staying, as shown in the previous paragraph, in the same component of 
T split by M'. Note that P' indeed uses only vertices of T(. This completes the induction step. ■ 



3.2 Approximation ratio analysis 

For a collection A of directed stars S(u i: r^), define w(A) = J2s( Ul ri )eA Ti > tne tota ^ P ower use d by 
the stars in A. 

Lemma 4 Let B be a an arbitrary collection of stars, and T be an arbitrary spanning tree. 
There exist non-negative coefficients x s (over the collection of all possible stars S(u,r)) such that 
Es x sfB(S)>c(T)-f(B)andZsXsp(S)<(l/2)opt 

Proof. We assume that c(T) — f(B) > 0, or else x s = for all S will do. Assign x s = (1/2) for 
every star of OPT, andrr^ = otherwise. Therefore J2 S xsp(S) = (1/2) opt. 

Recall that Q(B) = \J S ( Ui ri )et3Q( u ^ r i) an( ^ ^ et e e T \ Q(B). If we remove e from T, we 
create two subtrees T u and T v , where u and v are the endpoints of e. OPT, being strongly connected, 
has at least one star S u with the center in V(T U ) and one of its other vertices in V(T V ). Then e e 
Q[B U {S u }). Similarly, OPT has at least one star S v with the center in V(T V ) and one of its other 
vertices in V(T U ). Then e G Q(B U {S^}). Note that S v ^ S u (centers in disjoint vertex sets). 
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We have: 



s 



S e£Q{BU{S})\Q(B) 



E C M E 



eeT S | e6Q(BU{S})\Q(8) 



Ce XS 



e£T\Q(B) S | eGQ(S) 



> 



^ c e = c(T) - /(£) 



e£T\Q(B) 



where the inequation is given by the two distinct stars S u and S v described above for every edge 
e E T \ Q(B). This completes the proof. ■ 

The example from Figure [T] shows that a constant better than 1/2 is not possible above: with B 
being empty, we have c(T) circa n, and for each star S of optimum, p(S) is circa 1 and f(S) circa 2. 

Now we need the following lemma, whose proof is obtained from Robins-Zelikovsky as presented 
in 021 by changing what quantities represent and some parameters, together with a "fractional cover" 
idea from the journal submission version of [j3j . 

Lemma 5 Assuming that for the minimum spanning tree T, and for any collection of stars B, there 
exist non-negative coefficients (xs) such that J2s x sfB(S) > c(T) — f{B) and xsp(S) < aopt, 
where opt is the power of the optimum solution and a < 1, the algorithms' output has power at most 
flopt where (3 = 1 + a + a ln(l/a). 

Proof. First, if c(T) < aopt, then before any improvement we have a solution of cost at most 2a opt 
and 2a < (3. Thus in the following we assume opt > c(T) > aopt (the first inequality is due to T 
being a minimum spanning tree). 

Note that at the end of the algorithm, M contains exactly one of the two antiparallel arcs for each 
edge of T. Then, for the final collection of stars A, the output H satisfies 



which holds for every vertex u (recall that E(H) = [j SeA E(S) U M). 

Let Si, S 2 , ■ ■ ■ , S q be the stars picked by our algorithm and let Ai, for 1 < i < q be the collection 
of the first i stars; also let for convenience Ao be the empty collection. For 1 < i < q, let pi = p(Si), 
and let f { = f^iSi). Note then that for all i, since /^(Si) = f(Ai) - /(A-i), we have 



If fi = 0, then pi = and Equation [2]below holds. Otherwise, the greedy choice of the algorithm 
and the assumptions of the theorem for B = Ai-i give: 



p{H) < c(T) + w{A) 



(1) 



as it follows by summation over u G V from 




Pi < ft 



aopt 



(2) 



cm - e;=i fi 
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Define the function g : [0..c(T)] — > [0..1] by g(x) = aopt/(c(T) — x) for x < c(T) — aopt, and 
g(x) = 1 for x > c(T) — aopt. Then from Equation |2] and Fact |2] (that p^ < fi), we obtain: 

fi 

Pi < g(x)dx. 




c(T) - aopt 



Figure 6: The function g(x) is given by the solid curve. YH=iP% * s tne shaded area, in rectangles 
of width fi and height a ^Li ■ In this particular picture, aopt = (2/3) c(T) and therefore the 

integral is circa 0.94c(T). 
Therefore (see Figure©: 



* MT) j-c(T)—aopt ^ f c(T) 

2_Pi < / g{x)dx= I — — _ ^ cto + / _ Ida; 



i=i 



c(T)-»op< 



c(r)-aopt 

= (-aopt) ln(c(T) — x) + (c(T) - (c(T) - a op*)) 

1 + In 

Using this and c(T) < opt and Equation[T](recall that w(A) = Y^l=i P«)> we obtain that the power 
of the output is at most 

c(T) + aopt ( 1 + ln-S- J < opt (1 + a + cdn(l/a)) 
\ aopt / 

finishing the proof. ■ 

Based on Lemmas [3] and @J Theorem \T\ follows immediately from the fact that a — 1/2 makes 
/3 < 1.85. Note also that a < 1 implies /3 < 2, which follows from a(l + ln(l/a)) < 1, which is 
equivalent to ln(l/a) < l/a — 1, a fact that holds for all a < 1. HI had a = 7/8. 
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4 Linear Programming Relaxation 



While not improving the approximation ratio of Greedy, this may allow for further LP-based algo- 
rithm. The following natural Integer Linear Program is called IP2 as in [[5]|. We adapted the notation, 
and have variables ys for every arc S = S(u, r). The idea is that ys being 1 "represents" that S is star 
of the optimum solution. We say that star S = S(u, r) e 6~(X), for X C V, ^ X V iff u £ X 
and V(S) n X ^ 0. 



minimize 




ys> 1 vici/,«M^ (3) 

Ses-(X) 

ys>0 WS (4) 

ys e Z V 5 (5) 

LP2 is the linear relaxation of IP2; that is, the linear program given by exactly the same constraints 
except the last one. LP2 has exponentially many "cut" constraints, but when we replace them by 
"flows", we have 0(mn 2 ) non-zero entries in the matrix, as shown in [0. Thus it can be solved in 
polynomial-time . 

Let opt* be the optimum of the linear program LP2 for a given instance. Then clearly opt* < opt 
and [5J proves that opt < 2opt*. We do better here. First, Lemma |4]has a straightforward adaptation: 

Lemma 6 Let B be a an arbitrary collection of stars, and T be an arbitrary spanning tree. 
There exist non-negative coefficients xs (over the collection of all possible stars S(u,r)) such that 
E S xsfB(S)>c(T)-f(B)and>~:sXsp(S)<(l/2)opt*. 

Proof. We assume that c(T) -f(B) > 0, or else x s = for all S will do. Let y s (for all S = S(u, r)) 
be a an optimum solution of LP2. Assigns = (l/2)ys for all S; therefore J2 S xsp(S) = (1/2) opt*. 

Recall that Q(B) = \J S ^ U . r .) eB Q(ui, r/) and let e e T \ Q(B). If we remove e from T, we create 
two subtrees T u and T v , where u and v are the endpoints of e. Constraint |3] gives 

Se8-(V(T U )) 

and 

Ses-(V(T V )) 

Note that 5~ (V(T U )) and 5 _ (V(T U )) are disjoint sets since a star in the first set has its center in V(T V ), 
while a star in the second set has its center in V(T U ). Moreover, any star S G 5~(V(T U )) Li5~(V(T v )) 
has eeQ(S). Thus 

S | egQ(5) 
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From now on we coppied from the previous proof: 

S S e£Q{BU{S})\Q(B) 

= 5Z c ( e ) Xs 

eST 5 | e&Q(BU{S})\Q(B) 

e£T\Q(B) S | eeQ(S) 

> £ c e = c(T)-f(B), 

e£T\Q(B) 

where the inequation is from ©. This completes the proof. ■ 

The example from Figure [T]again shows that a constant smaller than 1/2 is not possible in Lemma 
with B being empty, we that for each star S with f(S) > e, p(S) > (1/2 - e)f(S). Thus, if 
EsXsf(S) > c(T), then J2 s x s p(S) > (1/2 -e)c(T) > (1/2 - e) (1 - 7 ) opt > (1/2 - e - j)opt*, 
where we used that the example can be made to have opt(l — 7) < c(T), for any 7 > 0. As 7 and e 
can be made arbitraly small, we can see that indeed a constant smaller than 1 /2 is not possible. 

It is also proven in Section II of [|5l that c(T) < opt*; then Lemma |5] goes through with opt* 
instead of opt. Therefore we conclude that the output of Greedy is at most 1.85 opt*. 

5 Conclusions 

This work greatly simplifies and at the same time improves our earlier work [4]. Instead of Greedy 
we could have used the the iterative randomized rounding of Byrka et. al for Steiner Tree, with 
the same approximation ratio. However, we do not see further improvements coming from using their 
full range of techniques, since we do not see the equivalent of the concept of "loss" used explicitly by 
[fT9l and implicitly by 0. 
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